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We analyze the equilibrium transport properties of underscreened Kondo effect in the case of a 
two-level quantum dot coupled to ferromagnetic leads. Using the numerical renormalization group 
(NRG) method, we have determined the gate voltage dependence of the dot's spin and level- resolved 
spectral functions. We have shown that the polarization of the dot is very susceptible to spin 
imbalance in the leads and changes sign in the middle of the S — 1 Coulomb valley. Furthermore, 
we have also found that by fine-tuning an external magnetic field one can compensate for the 
presence of ferromagnetic leads and restore the Kondo effect in the case oi S = \ Coulomb valley. 
However, the underscreened Kondo effect cannot be fully recovered due to its extreme sensitivity 
with respect to the magnetic field. 

PACS numbers: 72.25.-b, 73.63.Kv, 85.75.-d, 72.15.Qm 



I. INTRODUCTION 

Ever since the observation of the resistivity anomaly 
in normal metals at low temperatures/ the Kondo prob- 
lem^ has been investigated perpetually both experimen- 
tally and theoretically — As a result of that intensive 
research, in many respects, the Kondo effect is well un- 
derstood now. Very recently the interplay between the 
Kondo effect and other many-body phenomena, such as 
superconductivity or ferromagnetism, has attracted a lot 
of interest. That attention was mainly motivated by the 
recent advances in nanofabrication which have opened 
the possibility to attach superconducting^ii or ferromag- 
netioSi — leads to molecules or semiconducting quantum 
dots. In the following we will focus our attention on the 
case of ferromagnetic leads coupled to a quantum dot 
exhibiting the Kondo effect. 

From the theoretical side, a consensus was found that 
single level quantum dots (i.e. relatively small dots with 
level spacing 5 larger than the hybridization V) attached 
to ferromagnetic leads exhibit finite spin asymmetry as 
a result of the spin imbalance in the leads, which gives 
rise to a splitting and suppression of the Kondo reso- 
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This spin asymmetry can, however, be com- 
pensated by means of an external magnetic field or by 
tuning the gate voltage properlyJ^iiS In that way the 
strong coupling Kondo fixed point can be reached, with 
a fully developed Kondo resonance and a somewhat re- 
duced Kondo temperature. This expectation has recently 
been confirmed experimentally^"— 

While in many cases the quantum dots can success- 
fully be modelled by only a single local level, there are 
few exceptions when the finite level spacing affects the 
low energy physics of such a device considerably. An iso- 
lated quantum dot with even number of electrons can 



have spin triplet S = 1 ground state, if the level spac- 
ing between the two single particle levels closest to the 
Fermi energy is anomalously small, i.e. smaller than the 
exchange coupling between the two electrons situated on 
those levels. If such a device is attached to only one lead 
then the system will exhibit the so-called underscreened 
Kondo effect — A single lead can only screen half of 
the local spin below the Kondo temperature Tk while 
the residual spin-a-half object is left unscreened with a 
ferromagnetic Kondo coupling to the rest of the conduc- 
tion electrons. Such a ferromagnetic Kondo coupling is 
known to be irrelevant in renormalization group sense, 
as it scales to zero when the energy scale is lowered, re- 
sulting in a Fermi liquid plus a decoupled ^S* = ^ object 
at zero temperature. This behavior is reached in a non- 
analytic fashion resulting in singular terms in e.g. the 
single particle scattering rate. Such kind of behavior was 
recently termed as singular Fermi UquidS^^ 

In a real quantum dot subject to a transport experi- 
ment the situation is slightly more complex as the dot 
has to be coupled to two leads to drive current through 
it. Even if the second lead is just a weakly coupled probe, 
its presence introduces a second energy scale, <C T^, 
at which the residual spin-a-half is screened by the other 
lead mode. The screening takes place in two stages and 
the separation between the energy scales Tk and can 
be tuned by the asymmetry between the coupling to the 
leads. Since the Kondo temperature is exponentially sen- 
sitive to the coupling, a realistic value of asymmetry is 
enough to separate the two stages completely. In such a 
case, if the other experimentally relevant energy scales, 
such as temperature T or magnetic field 5, lie in between 
the two Kondo temperatures <^ B^T < Tk, then 
the system is well described by the underscreened Kondo 
model. In fact, very recently, underscreened Kondo effect 



was observed experimentally in molecular quantum dots 
coupled to nonmagnetic leads. 

In the present paper we focus our interest on the inter- 
play of underscreened Kondo effect and itinerant electron 
ferromagnetism in the leads. To that end we consider a 
two level quantum dot coupled to a single reservoir of 
conduction electrons exhibiting spin imbalance. We as- 
sume that the second Kondo temperature is always much 
smaller than the experimental temperature, therefore we 
consider the second lead as a weakly coupled probe only. 

II. THEORETICAL DESCRIPTION 




FIG. 1: (color online) Schematic of a two- level quantum dot 
coupled to a ferromagnetic reservoir with spin polarization 
p. The dot level energies are denoted by Si and £2, S is the 
level spacing, while U' and J describe the inter-level Coulomb 
correlations and spin exchange interaction. The jih dot level 
is coupled to ferromagnetic electrode with strength Fja. 



The considered system consists of a two-level quantum 
dot coupled to a ferromagnetic reservoir, see Fig. [TJ and 
its Hamiltonian is given by 



H = Hfm + ^QD + ^tun, 



(1) 



where the first term describes the noninteracting 
itinerant electrons in ferromagnetic lead, H^m = 
^ka^kacl^Cka, whcrc 4^ is the electron creation op- 
erator with wave number k and spin a, while Ska is the 
corresponding energy. The second part of the Hamilto- 
nian describes a two-level quantum dot and is given by 

jcr j 

-^U' J2 nian2a' + JS^ + BS', (2) 

where rija- = djcr^ja creates a spin-a electron 

in the jth level {j = 1,2), Sj denotes the corresponding 
energy of an electron in the dot. Here, ei = e — S/2 
and £2 = s -\- with S = €2 — Si being the level spac- 
ing. The on-level (inter-level) Coulomb correlations are 
denoted hy U {U')^ respectively. J is the ferromagnetic 
exchange coupling (J < 0) with S = 6'i + 6*2, where 
Sj = \ ^cr(T' ^\cr^crcr'djcr' Is the splu opcrator for elec- 
trons in the dot level j and <j denotes the vector of Pauli 
spin matrices. The last term of i^QD describes the Zee- 
man splitting with B {gjiB = 1) being the external mag- 
netic field applied along the zth direction. 
Finally, the tunnel Hamiltonian is given by 



kja 



-ka^J 



djcr 



(3) 



where tja- describes the spin-dependent hopping matrix 
elements between the jth dot level and ferromagnetic 
lead. The strength of the coupling between the dot and 
lead can be expressed as Tja- = 7rp|tjcrp, where p is the 
density of states (DOS) in the lead. In the following, 
we assume a flat band of width 2D and use D = 1 
as energy unit, as not stated otherwise. Note that by 
assuming constant DOS the whole spin-dependence has 
been shifted into the coupling constants. While this as- 
sumption may simplify the calculations, it does not affect 



the low-energy physics we are interested inJ^iii To pa- 
rameterize the couplings, it is convenient to introduce 
the spin polarization of ferromagnetic lead, p, defined as 
p = (r^ — r|)/(r^ + r^j^ where we have assumed that 
each dot level is coupled with the same strength to the 
lead, Tja- = F^r. Then, the coupling for the spin-^ (spin- 
D electrons can be written as T^^^) = (1 ± p)r, where 
r = (r^ + r|)/2. For p = 0, the couplings do not de- 
pend on spin and the system behaves as if coupled to 
nonmagnetic lead. However, for finite spin polarization, 
p 7^ 0, the couplings are spin-dependent, giving rise to an 
effective exchange field which may spin-split the levels in 
the dot, suppressing the Kondo resonance. Such behav- 
ior has been extensively studied theoretically in the case 
of single- level quantum dots^— i^^— where the usual 
/S = I Kondo effect develops, while the effect of ferro- 
magnetism on the other types of Kondo effect remains to 
a large extent unexplored. 

In this paper we thus consider the zero-temperature 
equilibrium transport properties of two-level quantum 
dots coupled to ferromagnetic lead in the case of under- 
screened Kondo effect. In order to perform this analysis 
in most accurate and exact way, we employ the numeri- 
cal renormalization group (NRG) method^^i^ The NRG 
consists in a logarithmic discretization of the conduction 
band and mapping of the system onto a semi-infinite 
chain with the quantum dot sitting at the end of the 
chain. By diagonalizing the Hamiltonian at consecutive 
sites of the chain and storing the eigenvalues and eigen- 
vectors of the system, one can accurately calculate the 
static and dynamic quantities of the system. In partic- 
ular, numerical results presented here were obtained us- 
ing the flexible density-matrix numerical renormalization 
group (DM-NRG) code, which allows to use arbitrary 
number of Abelian and non-Abelian symmetries.^ — In 
fact, exploiting symmetries as much as possible is crucial 
in obtaining highly accurate data. Because in the case 
of ferromagnetic leads the full spin rotational invariance 
is generally broken, in the case of finite spin polarization 
we have used the Abelian symmetries for the total charge 
and spin, while in the case of p = we have exploited 
the full spin SU (2) symmetry. 

Using the NRG, we have calculated the expecta- 
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TABLE I: The respective eigenstates \Q^S^) and eigenvalues 
Eq^s^ of the decoupled quantum dot Hamiltonian. Here, Q — 
(X)jcT ^3"^ ~ 2) charge and zth component of 

the spin in the dot, while IX1X2) denotes the local states with 
Xj = 0, t, d for zero, spin- up, spin- down and two electrons 
in the level j. 



Tl 






1 
1 


-2,0) = |00) 

— 1 —-\ — 1 1 n\ 


n 
u 


2 


^ _ S _^ 3J_ _ B_ 


Q 


-i,-|> = |o;> 


£ -\- 1 -\- £l — R 


4 


_ 1 1\ _ 1 A r)\ 

_ 1 1\ _ IQ ^\ 

-■-5 2 / ~ 1^ 1 / 


p+ 1 + M + 1 
^ 1 2 ' 4 '2 


5 


6 


0,0) = |dO) 


2£ - ^ + [/ 


7 


0,0) = |0d) 

o,o) = ^(in>- 


2£ + ^ + [/ 


8 


;t)) ^e + U' 


9 


0,-1) = |U> 

0,0) = ^(in> + 


2e + U' + 2J - B 


10 


;t)) 2£ + [/^ + 2j 


11 


0. 1) = 1 tt) 

1, -i> = |d;> 
i,-i) = Ud) 


2e + U' + 2J + B 


12 


3£ - 1 + ^ + [/ + 2[/^ - f 


13 


3£ + 1 + ^ + [/ + 2[/^ - f 


14 


l,|> = |dt) 

i,i> = ltd) 


3£ - 1 + ^ + [/ + 2[/^ + f 


15 


3£ + 1 + ^ + [/ + 2[/^ + f 


16 


2,0) = |dd) 


4.e + 2U + 4JJ' 



tion value of the dot's spin as well as the spectral 



function of the dot, Ajji(j{(jj) 



G^,j^{uj)\^ where Gj^/^(c(j) denotes the Fourier trans- 
form of the dot retarded Green's function, G^j,^{t) = 
—iQ{t){{dj^{t)^d}-,^({))}). The diagonal elements of the 
spectral function Aj^{uj) = Ajj^iuo) are related to the 
spin-resolved density of states, whereas the off-diagonal 
elements Ajj>crioo) may be associated with processes of in- 
jecting and removing an electron at different sites. The 
symmetrized spin-resolved spectral function of the dot 
can be found from 




-0.8 -0.4 

e/D 

FIG. 2: (color online) The expectation value of the dot spin 
operator 5*^ in presence of nonmagnetic (a) and the expecta- 
tion value of dot spin operator for ferromagnetic (b) lead 
as a function of the average level position e — [si -\- S2)/2. 
The solid line corresponds to F = 0.04, while the dashed line 
corresponds to F = 0.015 (in units oi D = 1). The expecta- 
tion values separately for the two dot levels are also shown. 
The parameters are 6 0.1, J -0.2, U = U' = 0.4, B = 0, 
and in the ferromagnetic case p = 0.4. The shadowed regions 
indicate Coulomb valleys with different electron number in 
the dot, where N — ^A^rtja- 



A,H = ^A,-,-vH. (4) 

33' 

On the other hand, the normalized full spectral function 
is given by ii^^TcrAcriuo) and is directly related to the 
conductivity of the dot. 



III. NUMERICAL RESULTS 

Before presenting numerical results, it is instructive 
to analyze the eigen-spectrum of the decoupled quan- 
tum dot Hamiltonian, Eq. (|2]). The respective eigenstates 
\Q^S^) and eigenvalues Eq^s^ of ^QD are listed in Table 
I, where Q = (Xljcr ^jcr ~ 2) denotes the charge in the dot 
and is the 2:th component of the dot spin. By tuning 
e the ground state of the dot changes. In particular, in 
the absence of magnetic field and for ferromagnetic ex- 
change coupling (J < 0), when e > £-2,-1 = f ~ 
the dot is empty, for e = £-1,0 = ~f ~ ~ X' 
doublet \Q = -1,S = |) and triplet |Q = 0,6' = 1) 



states become degenerate, where now S is the full dot's 
spin, whereas for e = £0,1 = ^ — U' — U-\-^, the states 
IQ = 0, 6 = 1) and |Q = 1, 6 = ^) are degenerate. Con- 
sequently, for 6-2,-1 > s > 5-1,05 ^he ground state of the 
decoupled dot is the doublet and the system will exhibit 
usual S = ^ Kondo effect, whereas for £_i,o > £ > £0,1 
the dot is in the triplet state and the system will show 
an underscreened Kondo effect. In calculations we have 
taken the following parameters: S = 0.1, J = —0.2, 
U = U' = 0.4, and F = 0.04 (in units of D = 1). One 
then finds, £-2,-1 = 0.2, £_i,o = —0.2 and £0,1 = 1- 

We have calculated the expectation value of the dot's 
spin operators S'^ and in case of nonmagnetic and 
ferromagnetic leads, respectively. The results are sum- 
marized in Fig. [21 As the level position is lowered by 
sweeping the gate voltage, the dot is tuned through the 
Coulomb blockade valleys with electron number = 
1, 2, 3 and spin states 6 = 1/2, 1,1/2, respectively. When 
the dot is attached to a ferromagnetic lead, the spin 
asymmetry in the hybridization results in the spin split- 
ting of the dot level, thus leading to a spin polarized 



4 



ground state. Precisely in the middle of the N = 2 
Coulomb blockade valley, i.e. for e = (^-1,0 + ^o,i)/2 = 
— U^ — ^ = —0.6, the spin polarization changes sign, as 
to the right (left) hand side from that point the electron 
(hole) like virtual processes dominate in the renormaliza- 
tion of the level position. 

Note that the abrupt jump of the dot polarization is 
due to the fact that the underscreened Kondo model is 
extremely susceptible to even a tiny magnetic field. The 
reason is that the ground state of the system consists 
of a Fermi liquid and a decoupled residual spin S = ^ 
and that residual spin at zero temperature can be polar- 
ized by any infinitesimal magnetic field. In the case of 
a weakly coupled second lead there is a second stage of 
the Kondo screening at when the remaining spin de- 
gree of freedom is quenched by the weakly coupled lead's 
electronsj^ In that case the sudden step in the polar- 
ization turns into a very steep crossover with a width of 
~ max(T^, T), where T is the experimental temperature. 

Since the dot spin polarization is rather difficult 
to detect experimentally, we have computed the level- 
resolved single particle spectral density Ajj'c^{(j) = 
-^lm[Gf^,^{uj)^Gp^^iuj)], where Gf^.^{uj) denotes 
the corresponding retarded Green's function. This quan- 
tity contains the information about the transport prop- 
erties of the dot. Given that the dot is strongly cou- 
pled to one of the leads, the conductivity through 
the setup at voltage bias eV is given by dl/dV ~ 

The results for the level-resolved normalized spectral 
functions -kTA^J^ = nJ^a'^aApl in the case of non- 
magnetic lead {p = 0) are shown in Fig. [3l The spectral 
functions are plotted as a function of energy uj and aver- 
age level position e. To resolved the low-energy behavior 
of spectral functions logarithmic scale for u is used. Fur- 
thermore, only the Coulomb blockade valleys with N = 1 
and N = 2 electrons in the dot are shown. The behav- 
ior of spectral functions for blockade valley with N = 3 
electrons is similar to that with a single electron due to 
the particle-hole symmetry. 

First of all, we note that the two Coulomb blockade 
regimes with = 1 and N = 2 are characterized by fun- 
damentally different fixed points. In the N = 1 regime 
(—0.2 < e < 0.2) the effective low-energy model is a 
Fermi liquid with fully screened spin S = ^ {S = ^ 
Kondo model). On the other hand, in the two-electron 
regime (— 1 < £ < —0.2) the effective model is a singu- 
lar Fermi liquid^ i.e. a normal Fermi liquid plus a free 
spin- 1 with weak residual ferromagnetic coupling to the 
conduction channel. These two distinct ground states 
of the system are separated by a quantum phase transi- 
tional^ which occurs when sweeping the gate voltage 
through the boundary of = 1 (A^ = 3) and AT = 2 
Coulomb valleys. The quantum phase transition is of 
Kosterlitz-Thouless type^ with exponentially decreas- 
ing Kondo temperature when approaching the transition 
from the normal Fermi liquid side. The different be- 
havior associated with two distinct fixed points can be 
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FIG. 3: (color online) The density plots of the normalized 
level- resolved spectral function TrFA^y,^ = TrJ^a^^^^^ft (^)" 

(c) and the full spectral function ttFA^^^ = E,,- '^^^^jj'^ 

(d) in the case when the dot is attached to a nonmagnetic 
electrode (p = 0). The parameters are the same as in Fig. [2] 
with r = 0.04. Note the logarithmic scale on the frequency 



observed in the dependence of the level-resolved spec- 
tral functions. In the N = 1 Coulomb valley the diag- 
onal elements of the spectral function exceed the uni- 
tary value of ttF, while the off-diagonal elements become 
negative and their absolute value is also larger than ttF. 
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FIG. 4: (color online) The density plots of the spin-dependent 
level- resolved spectral function nVaA-^y,^'^ (a)-(f), the spin- 
dependent full spectral function vrr^Ar^'^ = E^/ ^^^ij^a ^ 
(g) and (h), and the full spectral function 7vJ2a Ta^S^^'^ (i) 
in the case of ferromagnetic lead with p = 0.4. The param- 
eters are the same as in Fig. [2] with F = 0.04. Note the 
logarithmic scale on the frequency axis. 



On the other hand, in the N = 2 Coulomb blockade 



valley, the diagonal and oflF-diagonal elements of A^^^ 
are rather positive and always smaller than ttF. In con- 
sequence, the full spectral function is properly normal- 
ized = + AP=^ + 2Af J° = 2^F for = 0, 
see Fig. [3fd). Furthermore, the full spectral function 
clearly exhibits the resonance corresponding to the un- 
derscreened Kondo effect in the Coulomb blockade valley 
N = 2{e< -0.2), while for = 1 > -0.2) the S = ^ 
Kondo effect develops. 

The spin and level-resolved normalized spectral func- 
tions 7tTA^~,^'^ in the case of ferromagnetic lead with 
p = OA are displayed in Fig. SI Now the spectral func- 
tion is different for each spin component due to the spin- 
dependence of the coupling parameters. Furthermore, 
when the lead is ferromagnetic, both Kondo resonances 
become completely suppressed. This is directly associ- 
ated with spin-dependent level renormalization, which 
lifts the spin degeneracy in the dot, suppressing the spin- 
flip cotunneling processes driving the Kondo effect. In 
other words, the ferromagnetic lead exerts an effective 
exchange field on the dot, which destroys the Kondo ef- 
fect. There is only a very narrow resonance at e = —0.6 
when the spin splitting of the dot levels generated by 
the spin imbalance of the conduction electrons vanishes, 
which happens exactly in the middle of the triplet valley, 
see also Fig. EJb). 

As mentioned above, the suppression of the Kondo 
peaks is associated with the exchange field which devel- 
ops in the presence of spin-dependent couplings. How- 
ever, as shown in the case of single level quantum 
dots^ii^ one may try to compensate for the presence of 
the exchange field by applying properly tuned external 
magnetic field. In Fig. [5] we show the spin-resolved spec- 
tral functions calculated in the regime of S = ^ Kondo 
effect for £ = for different values of magnetic field B. It 
can be clearly seen that by tuning the magnetic field it is 
possible to fully restore the Kondo resonance at = 0. 
This happens when B = = 0.009215, see Fig. El with 
Be being the compensating field, i.e. a field at which the 
previously-split dot level becomes degenerate again. 

The situation, however, becomes more complicated for 
e = —0.4, which corresponds to the Coulomb valley with 
N = 2 and 5 = 1, when the dot is described by the un- 
derscreened Kondo model. The respective spin-resolved 
spectral functions are shown in Fig. [6] using the loga- 
rithmic scale to emphasize the distinct dependence on 
magnetic field. Because the dot is coupled only to one 
conduction channel, only a half of the dot's spin can be 
screened by conduction electrons and the ground state 
consists of a Fermi liquid and a decoupled spin S = \. 
Consequently, at T = an infinitesimally small mag- 
netic field is enough to polarize the residual spin-a-half 
in the dot. This leads to an extremely high sensitivity 
of transport properties with respect to magnetic field, as 
can be seen in Fig. [5l In consequence, in order to com- 
pensate for the presence of exchange field in the case of 
underscreened Kondo model, one needs to perform a fine- 
tuning of magnetic field. Furthermore, it turns out that 
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FIG. 5: (color online) The spin-up (a), spin-down (b) and the 
full (c) spectral function in the S — ^ Kondo regime for e = 
in the presence of external magnetic field B, as indicated in 
the figure. The parameters are the same as in Fig. [2] with 
r = 0.04 and p = 0.4. The inset in (c) shows the zoom- 
out full spectral function. The Kondo effect is restored when 
B = 0.009215. 



FIG. 6: (color online) The spin- up (a), spin-down (b) and the 
full (c) spectral function in the underscreened Kondo regime 
for s = —0.4 in the presence of external magnetic field B, as 
indicated in the figure. The parameters are the same as in 
Fig. [2] with r = 0.04 and p = 0.4. The inset in (c) displays 
the dependence of the spectral function associated with the 
Kondo peak. 



although it is possible to restore the resonance peak at 
cj = 0, the full restoration of the underscreened Kondo 
effect is not possible, as the height of the peak is below 
the unitary limit. 



IV. CONCLUSIONS 

In this paper we have analyzed the equilibrium trans- 
port properties of a two-level quantum dot asymmetri- 
cally coupled to ferromagnetic leads. We have shown 
that by tuning the position of the dot levels, the ground 
state of the system changes from a Fermi liquid in the 
S = ^ Coulomb valley into a Fermi liquid plus residual 
spin-a-half in the S = 1 Coulomb valley, the latter being 
the example of underscreened Kondo effect. The bound- 
ary between these two regime is a quantum phase tran- 



sition. For finite spin polarization of the leads, p 7^ 0, 
the Kondo phenomenon becomes suppressed due to an 
effective exchange field originating from the presence of 
ferromagnetic leads, and so is the critical behavior. 

In the transport regime where the system exhibits un- 
derscreened Kondo effect and for p 7^ 0, we have found 
that the polarization of the dot abruptly changes sign 
in the middle of the S = 1 Coulomb blockade region. 
This is associated with virtual tunneling processes (either 
electron-like or hole-like) that give the dominant contri- 
bution to the renormalization of the dot levels. 

Furthermore, we have also analyzed the effect of an 
external magnetic field B applied to the dot. It has been 
shown that an appropriately tuned B can restore the 
Kondo effect in the case of the S = ^ Coulomb valley. 
The underscreened Kondo effect, on the other hand, ex- 
hibits an extremely high sensitivity towards even a tiny 
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change in magnetic field due to its particular ground 
state. Consequently, an infinitesimally small magnetic 
field is enough to polarize the residual spin-a-half in the 
dot, which definitely hinders a full restoration of the un- 
der screened Kondo effect by tuning the magnetic field. 
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